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Abstract 



We introduce one matrix model coupled to multi- flavor vectors. The two- flavor 
vector model is demonstrated to reproduce the two-point correlation numbers of 
boundary primary fields of two dimensional (2, 2p + 1) minimal Liouville gravity on 
disk, generalizing the loop operator (resolvent) description. The model can properly 
describe non-trivial boundary conditions for the matter Cardy state as well as for 
the Liouville field. From this we propose that the n-flavor vector model will be 
suited for producing the boundary correlation numbers with n different boundary 
conditions on disk. 
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1 Introduction 

There have been developed two independent approaches to (Euchdean) two dimensional 
quantum gravity since the middle of 80's: Liouville gravity (LG) P-S] versus matrix 
models (MM) [SHlOj. Two approaches are checked in a number of particular models 
[TTVfTS] . One may refer to more references in reviews in [T^[T5] . 

Initiated by Al. Zamolodchikov's direct computation of partition function on fluctu- 
ating sphere at the simplest integral point (Liouville coupling set to b"^ = |) [16], Belavin 
and Zamolodchikov succeeded to confirm the correlation numbers of minimal gravity by 
providing so-called resonance transformation formula between conformal (Liouville) and 
KdV frame [I7j. 

Correlation numbers are given in terms of the integrated form of the local density (2- 
form operator) Ok{X), Ok = jj^ Ok{X) over the manifold Ai which accommodates both 
"matter" and the metric degrees of freedom localized aX X & M.. A generating function 
of the correlation numbers in Liouville gravity is given by ^lg({'^A:}) = (6^'='^'=*^'= ) . A 
similar generating function, which depends on the parameters tk-, can be introduced in 
Matrix models. It was conjectured in [18] that there exists a special choice of contact terms 
in MM or, equivalently, the special transformation tk = tk{{Xj}) such that the coincidence 
of the partition functions is ensured in the p-critical one-matrix model OMM(p) and the 
Minimal Liouville gravity MLG(2,2p+l) for random surfaces with spherical topology. This 
relation between the parameters tk and was obtained in [18] up to the linear terms. Its 
explicit form to all orders was conjectured in [17] and checked up to the 4-th order [TTlflQ] . 

Given the minimal gravity, it has been also shown [20] that the resonance transforma- 
tion formula also works for bulk correlation numbers in the presence of fluctuating disk. 
In particular, the disk partition function of OMM(p) with the boundary length I in the 
conformal frame is given as [T8|[20] 



where /i is the bulk cosmological constant. On the other hand, the Liouville gravity 
partition function Zd^^.^b) is related by the inverse Laplace transform of Zc){fi,l) 




(1.1) 




(1.2) 
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where the contour f goes along the imaginary axis to the right from all the singularities 
of the integrand. Using the Laplace transform one can show that the resolvent in the 
continuum limit on the spherical topology 



is proportional to cosh(7rs/6) when the boundary parameter is set to z = cosh(7r6s). 
Mo provides a scale parameter in the matrix model (details are given in appendix A). For 



In Liouville field theory (LFT), b corresponds to the Liouville coupling constant relat- 
ing the background charge Q = b + l/b and z to the boundary cosmological constant fiB = 
y^7'7shi(^'^^cosh(7r6s). MLG(2,2p+l) corresponds to the case b^ = 2/(2p + l). Boundary 
operator of LFT, [e^'^{X)'j '^^ with X E dAi has conformal dimension A(/3) = l3{Q — /3) 
and is specified by the boundary condition BC(si; S2) which e^"^ joins. The primary Liou- 
ville boundary operator in MLG(2,2p+l) is given as Bu = e^^^^'^u and l3u = b{l + £)/2 
with i = 1,2, ■ ■ ■ ,p. $mn represents the CFT matter in the Kac table with the conformal 
dimension A^n = C(mn{c(mn — 1) where a^n = {n — l)b/2 — (m — l)/(26) and q = 1/6 — 6. 
It should be noted that in LFT, the boundary correlation of e^'^'s is known in [2n[2l] . 
not to mention the boundary correlation for the A-series minimal model [22l[23]. How- 
ever, the corresponding result is not available in the matrix model side, even though few 
attempts can be found in RSOS model and 0(N) fluctuation models [25l[26], in loop gas 
model [271EBJ, and in (two) matrix model [291180]. 

In this paper, we propose boundary changing operator description in OMM(p) which 
reproduces boundary correlation numbers of primary fields in MLG(2,2p+l) on disk. This 
paper is organized as follows. In section [21 we introduce one matrix model which couples 
to certain number of vectors. In section [3l we demonstrate that 2-point correlation number 
of the integrated form of boundary changing operator, = ^■^'^'^^ B-I'p'^\x) 

are reproduced from the matrix model with two vectors. Here Si and S2 correspond 
to the boundary conditions of LFT and 1 = (1,1) and i = are the Cardy label 

of minimal boundary conditions. Starting from BC(si,l; S2,i) one can construct the 
general boundary condition, BC(si, m; S2, n), when 1 < m,n < p are allowed from fusion 
property. Section H] is devoted to summary and discussion. 




(1.3) 



p-critical theory, u{z) has the scale dimension Uq so that u{z) 



,P+l/2 



cosh(7rs/6). 
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2 One matrix model with vectors 



In order to describe 2 dimensional gravity with boundaries, we introduce one matrix model 
with vectors, 



J DMBv^'^^^Dv^^^ exp ^-^trr(M) " Z]^^"^^ " C^'^'^'Km) ■ v^^^^ , (2.1) 



where M is a x Hermitian matrix and V{M) is a polynomial of M which starts 
from a quadratic term . f and its hermitian conjugate v'^°''>^ are A^ dimensional 
vectors and ■ represents the contraction of the A^ dimensional indices, a and h label the 
"flavors" of the vectors. The number of flavors may depend on the number of different 
boundary conditions. C*-"'*-'(M) is a hermitian matrix whose diagonal component is given 
as a polynomial of M , 

CM(M) = ^fo^f^M^"-'^, (2.2) 



with 6q"^ normalizecjl 



n=0 

do) 



as (—1)^" and K"^ is the order of the polynomial, hn in the diagonal 
component behaves as the source to the boundary preserving operator. The off-diagonal 
component is also given as a polynomial in general and its coefficients are the source to 
the boundary changing operator. 

This model can generate Feynman diagrams which correspond to 2 dimensional dis- 
cretized surfaces. We associate a double line with a propagator of M and a single line 
with that of the vectors as usual. The rule of drawing the Feynman diagram is shown in 
Fig. [U A typical Feynman diagram with a boundary is found in Fig. |2l Note that the 
vector model fl2.ip is quadratic in the vectors and thus, the propagators of the vectors 
always form loops which are regarded as boundaries. 

The partition function Z in (12. ip can be expanded over various topologies including 
disks as following. We first integrate out the vectors and obtain 

= j DM exp (^-^trV(M) - TrlogC(M)^ = {e-^''°^^^^''^)o, (2.3) 



^ Note that by rescaling the vectors, we can always fix the coefiicient of the highest order term in 
M . Nevertheless, to compare the matrix result with the Liouville gravity one can choose a convenient 
normalization. More details are found in Sec. [3l 
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Figure 1: (a): Propagator for M. (b)(c): Examples of bulk interactions, (d): Propagator 
for the vectors. (e)(f): Examples of boundary interactions. 

where Tr stands for the trace over both flavor and matrix indices, e^" is the partition 
function without vectors 



= j DMexp (^-^tiV{M)^ (2.4) 

and describes manifolds without any boundary. (■ ■ ■ )o in (12.31) stands for a normalized 
expectation value with respect to (12. 4p . By taking the logarithm of (12.31) . we obtain 



oo 



Z = Zo + (e-T'-i^sCCAf)^^ ^Zo + J2 Zh, (2.5) 

h=l 

where (■ ■ ■ )c denotes the connected part of (■ ■ ■ )o. Z^, becomes the partition function with 
h > 1 holes (boundaries) 

Z, = i((-TrlogC(M))'^). = l (^[f f Tr e-'^w)'^^. (2.6) 

Noting that an insertion of a traced operator lowers the order of by 1 and thus, 
Zh has the extra factor A^"'^ compared with Zq, we can topologically expand Z^ = 
Yl'^=o ^h^^i where x = 2 — 2(7 — /i is the Euler characteristic for the manifold with g 
handles and h holes. In the following, we restrict ourselves to the case of disk with no 
handle [g = Q,h = 1) for simplicity. 




Figure 2: Typical Feynman diagram with a boundary. 

3 Correlation numbers from one matrix model 

Let us first consider the one-flavor case. Suppose C{M) = x — M. In this case, (12. 6 p with 
h = 1 reduces to the well-known simple partition function on disk, 

Z, = -(Tr log(a; - M)), = (|^" ^Tr g-'^^"^^)^^ . (3.1) 

The constant x couples to the length of boundary according to (11. 2p and is identified 
as the boundary cosmological constant, which is real. Note that differentiating Zi with 
respect to x gives the resolvent uj{x). 

Next, suppose the order of C(M) is 2 (iT = 2). 

C(M) = C2 + ciM + = (x+ - M){x- - M) . (3.2) 

Ci and C2 are real and x~^ + x~ = —c\ and x+x^ = C2. A certain correlation will be 

generated if one differentiates Z^ with respect to the source C2- 

^, , , / 1 \ w(x+) — u'(x~) 
0(x+, xr) = ( tr-— — — — ) = J ^ . (3.3) 



(x+ - M)[x- -M) I x+ - X 

When x^ are real, 0(xj[^,x]~) is not vanishing and becomes the two point correlation 
number of the boundary cosmological operator B\\ |30]. However, if x^ is allowed complex 
this is not the whole story since 0(x^,xj~) can vanislj^. In this case, dZ^jdc^ can be 
interpreted as one point correlation of i5i3 on diskjfl The vanishing condition is u;(x+) = 



^This requirement is reminiscent of the method of solving the resonance relation on sphere [171 . 
•^ (|3.3|) is obtained from p.7p when xi — >■ oo whose result is reproduced by fusing two Si2's. We thank 
Bourgine for this point. 
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u){x~), or ± = 2ibn with n integer. In the following, the notations will be used 
x^^") = Uocosh{-Kbs^) and w{x^^^) = Uq°^^^^ cosh{Tcs^ /b). 

One can generalize the above argument to the case when C{M) is the £-th order 
polynomial of M. C{M) is factorized in terms of ^ solutions, {x^^\ ■ ■ ■ satisfying 
C(x*'-^'') = 0. Vanishing requirement on the one-point correlation 0{x^^\ ■ ■ ■ ,x^^^) gives a 
constraint between any of x^^^^s so that u{x^^^) = u{x^^^) (1 < j,k < i); 

^^4^ = mod 1 . (3.4) 
2ib 

Next, two different boundary conditions need two flavors. We put C(M) as 



CiM)= Z... (3-5) 




where F{M) is a certain polynomial of M and c^^ is a M-independent constant and is 
the complex conjugate of c^^. Differentiating the partition function Zi with respect to 
the sources c^^ and c^^ one has two point correlation of boundary changing operators 



c'J =0 



When F{M) is linear in M, the correlation trivially reduces to that of -Bn's. To 
describe non-trivial ones, we may putF(M) = {xi^^-M){xi ^-M) and choose = S2± 
ib with S2 real. This choice ensures ui{x'^^) = ui{x'^^) = —oo{x2) with X2 = Mq cosh(7r6s2) 
and C{M) hermitian. As the result, (13. 6p becomes 

1 



0{xi; {x2}2) = { tr 



(xi - M){xi^^ - M){xi'^ - M) 
ul~^ cosh (^) cosh (^) 



2sinh ( ^M££±!^) sinh ( ^Mrf) ginh ( ^Mfi^) ginh ( ^MfiiizM 



(3.7) 



where Sp = Si + S2 and Sm = si — S2. This exactly agrees with Liouville two-point 
correlation of B12 with BC(si; S2) for the Liouville field. 

The generalization to other boundary operators is surprisingly simple: Put F[M) as 
the £-th order polynomial and choose (with S2 real) 

sl = S2 + tbj for j = -(£-l),-(£-3),--- ,£-3,£-l, (3.8) 
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Xl 




Xl 




Xl 




Figure 3: The boundary 2-point correlations for Bu, B12 and Bu. 
represents the factor of ^ 



The hne with x 



x-M' 



SO that (X'(x2"'^) = ui{x2^'''). Note that the number of allowed complex parameters for two 



.{i±2) 



point correlation is i as shown in Fig. [31 Explicit evaluation of fl3.6p is given as 

1 1 W , ,e uj{xi)-uj{xt'^) 



0{xi] {x2]i) = ( tr 



Xl - M Fe{M] X2] 



(^-l-2fc)N 



(3.9) 



where Fi{M; x) is the polynomials of M with degree i, Ylt^oi^^^"^'^''^ ~^)- This coincides 
with Liouville result of '^'^B^j up to BC-independent normalization [21]. (See (A. 16) of the 
first reference in [27] for the explicit formula for 6^ = 2/{2p + 1).) 

One can also specify the boundary condition for the matter field in (13. 9p . Noting 
the boundary structure constant c^^'^'^^^^ = 1 [23j, one concludes that Eg. (13. 9 p produces 
(BuBii) with BC(si, 1; S2, i) and Fi{M; x) generates i = (1, i) Cardy boundary condition. 
This conclusion is backed up by one-point correlation of Bn with BC(s, k), 



Oi{x}k) = ^trlog(Ffc(M;x)) = [k]guix) 



(3.10) 



where [x]q = {q^ — q^^)/{q — q^^) is the g— number with q = exp(i7r6^). From the field 
theoretic point of view, one can put the ratio 0{{x}k)/0{{x}i) as the ratio of the vacuum 
expectation value [22] of the matter part (1)^/(1)-,^ = [k]g since Liouville contribution 
cancels out in the ratio. The result is consistent with Eq. (I3.10p . 

Incidentally, we note that the matter operator is identified as = $i,2p+i-fc and 
BC(s; A;)=BC(s; 1, 2p+l—k) and Bi^k is related to Bi^2p+i-k through the Liouville reflection 
( "unitary" condition) [21]|27] . Since the matrix result Eq. (13. 9 p shows the same functional 
dependence as the Liouville result, it goes without saying that the same Liouville reflection 
holds for the matrix result [27j. Thus, even though the boundary primary operator Bi^k 
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is allowed as (1 < A; < 2p) and the order of the diagonal component of C(M) in Eq. (12. 2p 
can be restricted to 1 < K'^ < 2p for the p-critical model, the number of independent 
operators is further reduced by half by the Liouville reflection. 

From the above consideration, we assert that 0({xi}^; {x2}m) describes the two-point 
correlation of with BC(si,£; S2,m). One can support this idea using the fusion 

rule. Let us consider the four-point correlation with simple BC {i < m for definiteness) 

[B,,f'''^ [i5^,](^3;i) (s3;i) [i^^jC^^^-) (-4;m) [B.^J^''')^ . (3.11) 
If one uses the operator fusion rule 

{S4;m) TK? ]{si;l) ^ (si;l) r» ]{S2-/) _ rr, {si;m) \K? ](s2;^) 

[oim\ ^ [i^ie\ — ^\m-e+l\<k<m+e~l PlfeJ 

' [Ou\ <&) [Olm\ — ©|m-^+l|<fc<m+£-l ^ ' [tJlk] [6.iZ) 

one ends up with the two point correlation with general EC's. 

'{S2;i) (S4;m) . (3.I3) 

From the matrix point of view, the four point correlation can be conjectured of the 
form 

Two point correlation is obtained by contracting the 1-boundary parts, which expands 
l/Fi{M;x) in powers of M/xfl 

tr ( A'F— — i TTTTTT V ) ) (3-15) 



F,(M;x2)F^(M;x4 

with < ^ < £ - 1@ This shows that the number of allowed correlations is the same as 
that of Eq. (I3.13P and the range of k is i + m— 1 — 2g<k<l + 'm— 1. It is noted that 
(I3.15P with g can be generated by C(M) with M-dependent off-diagonal blocks. 



Fe{M;xi) G{M) 
G{My FUM;x2) 



ciM)^{ ^r;::/^ ..r::. i- (3.i6) 



When g = 0, Eq. (I3.15P is indeed the 2-point correlation of Bi^i+m-i with BC(si, S2, m) 
up to a normalization. To elaborate on this, let us consider the two-flavor matrix model 



M / x-expansion is equivalent to the small length expansion if one uses the Laplace transform Eq. (II. 2p . 
^ When g > i, Eq. p.lSp reduces to the summation of Eq. p.lSp 's with g < £ except £ — m = 1. 
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of the form fl3.16p with G{M) given by a M-independent constant: G{M) = c^^. For 
i = k = 2, one finds a relation after a short calculation, 

g^l^i =Oi{x^}2■,{x2}2)=2cos{nb'')0{xu{x2}3). (3.17) 

This relation connects BC(si, 2; S2, 2) and BC(si, 1; S2, 3) for the correlation of ^613. One 
can show a general recursion relation between the two-point correlations of Bi^e+m-i with 
different boundary conditions: 

[m - 1], O {{Xi}f, {X2}m) = O {{Xije+i, {X2}m-l) • (3.18) 

Thus, the matrix model predicts the obvious connection between BC{si,a; S2,b) with 
a + b = i + m for the correlation of Bi^i+rn-i- 

One can demonstrate that fl3.18p is consistent with the field theoretic result. For the 
matrix side, one can rescale the matrix Fk{M;x) so that it has the form, Nk Fk{M; x) 
where Nk is a constant. Then, the two point correlation is rescaled as 

^r^4^0({xi}fc;{x2W) ■ (3.19) 

For the field theoretic side, one considers the operator Ai^k+e e^^''''^'^'^ ^i,k+e, introducing 
the normalization constant Ai independent of the BC. The two point correlation of 
Bi^k+e with BC(si, k; S2, £ + 1) is given as 

^Ik+e dLi(3i,k+i, si, S2) dM{k; i + 1) . (3.20) 

Here, di is the Liouville correlation and d-M is the matter part correlation. Equating two 
results, Eq. (KT9^ and Eq. (K20^ one has 

O ({xi}fc; {X2}i+i) = Nk Ne+i Alk+e diWi^k+e, si, S2) dM{k; i + I) . (3.21) 

On the other hand, Eq. (13.211) should be compatible with Eq. fl3.18p . which leads to 
nontrivial relations between A^fc's. Nk Ng^i duik, £+1) = Nk+i Ni dM{k+l, i) . (Note that 
Liouville part cancels out because the same Liouville BC(si,S2) and the same operator 
Bi^k+e are used in both sides.) dM is given in terms of boundary structure constant [25] 

, iA _ (.X _ (£+l,M+l)(ll) (ooo] 

aMyi^,<- ^ -LJ — t^(^i^k+e),ii,k+£) \^/k — ^(i,fc+f),(i,fe+£) \^h+i ■ yo.zz) 
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This simplifies the relation as 

One can find a consistent solution of Nj^'s to Eq. fl3.23p which is crucial for the consistency 

of fl3.2ip . Considering the boundary structure constant is identified with the fusion matrix 

A, C 



[23] c 



{A,B,C) (K) 
LM 



B,K 



L, M 



and the fusion matrix is written in a factorized form 



when K = (11) [ST], F(i 



(Ik), (Ik) 
(Ij), (IJ) 



^ with Fi 



l/[i]g, one has 



fc+i 



AT. 



FfFi 



(3.24) 



FkFk+i 

Thus, Nk{k > 3) is determined completely from A^^i and N2. Finally, putting k = 1 
in Eq. fl3.2ip . one finds the field normalization A^^ in terms of NiNi. Once A^^'s are 
determined by the recursion relation f l3.24p . the identification (13.211) is consistent with 
the relation fl3.18p found in the matrix model. 



4 Summary and discussion 

We propose a generalized partition function of one-matrix model (12.31) to give boundary 
correlation numbers of primary fields of MLG(2,2p+l) on disk. We demonstrate explicitly 
that the two-fiavor vector model correctly reproduces the two-point correlation numbers 
of boundary primary operator Bi^k ^ k < 2p) in MLG(2,2p+l). 

A few comments follow. First, one can obtain the general boundary condition by 
contracting 1-boundary on the disk as given in Eq. (I3.13P and Eq. (I3.15p . One may 
equally contract i and m-boundaries in Eq. (13.111) . In this case one is left with 

(^^''^ [Buf'''^ ® ^'^''^ [Buf'''^ = ^''-'^^ [Bn]'^"'''^ (4.1) 

since (1, l)-boundary allows only (1,1) operator. This property is seen in the matrix side 
by contracting i and m-boundaries in Eq. (13.141) . Eq. (I3.15P is left with the two-point 
correlation of Fis only 

{f,{M;x,)F,{M;xs)) ) ^^'^^ 
10 



and its contribution is Bu with BC(si, 1; S2, 1). 

Second, the non-vanishing power oi M {g ^ 0) in Eq. 03.151) can be obtained using 
C{M) in Eq. ( 12. 2p when its off-diagonal component contains the M dependent term up 
to the power g which is less than any degree of the corresponding diagonal component. 
Therefore, the M-dependent off-diagonal term produces the two-point correlation of Bi^k 
of BC(si, a; S2, h) with a + h ^ k + 1. However, it is not clear yet, what kind of mechanism 
restricts the maximum power g. In addition, fusing rule is not simple for non- vanishing g 
as shown in Eq. fl3.15p . This suggests that one needs to find the more detailed description 
of the M-dependency of the off-diagonal term. 

Finally, is the proposal working for three point correlation with 3-flavor vectors? Let us 
consider 0{xi; X2; X3). This produces the 3-point correlation of iSn's in [30]. On the other 
hand, the Liouville boundary 3-point correlation C^'^^'^^g*^^ of e'^'^^'s with BC(si; S2; S3) 
satisfies the difference relation^ 



- sinh {f - si + z6 - 2^/33)} sinh {f (S3 -s^+ih- 21(5^)] C^^T^'''^ 
+ sinh {f (S2 + si+ih + 22/33)} sinh {f (S3 + si+ih + 2if32)] C^^XT'''^ 

= sinh {f (S2 + S3 + 4tb - 22/3i)} sinh(7rfoi)Cjf^+;'/^+^'''^^\ (4.3) 

0{xi] X2] Xs) obviously satisfies f l4.3p with f3i = (32 = (3^ = b. One can consider a more 
general form. For example, 0{xi; X2', {xs}2) satisfies the difference equation f l4.3p with 
(/3i,/32,/33) = (&, f, y) and corresponds to (-611-612-612) with BC(s2,2; S3, 1; Si,l). Like- 
wise, we expect that the n-point correlation number 0(xi; ■ ■ ■ ; x„) will produce the result 
(-B11 ■ ■ --Bii) with BC(si, 1; ■ ■ ■ ; s„, 1). Nevertheless, we need more careful check for cor- 
relations of other non-trivial operators of -Bi^^'s. In addition, it looks very plausible that 
the same idea can work for the boundary correlation of multi-matrix model. Further 
details will be reported elsewhere in the near future. 

Note added: After completion of this work authors of [33] let us know that the same 
model arises in string theory on partially resolved singularities. 



^We fix a typo in the original equation. 
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A Resolvent in the spherical limit 

We present an expectation value of the resolvent 

in the large- limit. In the double scaling limit (continuum limit), M ~ ai — where 
Oi, 02 are constants and Q = -^ + u{x) [71ll8y 32j. Thus, one may renormalize z — )■ 01 + 02^; 
and put 

z = uq coshi^nbs) (A. 2) 

In the large- A^ limit, u{x) is given by a maximal real solution to the string equation 

p-i 

P{u, {4}) = M^'+i + J2 tku^'"'' = 0, (A.3) 

fe=0 

with tp_i = X. The parameter —to is proportional to the bulk cosmological constant fi 
and the other parameters, {tk\l < k < p — 1}, describe the relevant deviations from 
the p-critical point in the KdV frame. They are related, through the so-called resonance 
transformation, to the perturbative coupling constants which couple to the operators 
in the minimal gravity. The resonance relation was solved in [1^ and P{u, {tk}) is given 
in terms of {Xk} as follows, 

P+i (P + 1)' f Lp+i{u/uo) - Lp_i(m/mo) 



Pin, {tkiXk)}) = ur [ ' - ' - J + 0(A.), (A.4) 

where uo is a solution to the string equation with tk = for A; > 1, namely, uq = y/—to ~ 
y/JI and Lp{x) is the Legendre polynomial. 

In the large- A^ limit, we can neglect the commutator [^, u{x)] = 0, 



2a2V7rZ 



to 

12 



Afe=0 



Jo 2a2Vnl Ju(x=tn) du 



Jo 2a2V 7r£ Ju{x=to) 



Afc=0 



roo 



2(2p-l)!!a2yo 7i 



If one uses 



f dxLr,{x)e = y — irp+i/2(0 

for the Macdonald function of a half-integer order i^p+i/2(a;), one obtains the singular 
part of the resolvent in the large- limit, 

uj{z) = auQ^^ cosh (^-6) 
where a — ^""^^2(2^4^1)!!^^^' independent of s and can be absorbed in the renormalization. 
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